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In this note one considers 2-diniensional identical particles with fractional statistics 
(anyons) [1]. This system seems to play an important role in the understanding of the 
Fractional Quantum Hall Effect [2]. 

The iV-anyon Hamiltonian (there is no additional mutual interaction) reads 

-, N AT , 

^^ = ^ E (pi -^Ajf + T. -^rn^'^r] (1) 

^^ 1=1 1=1 ^ 

where a is the statistical parameter (a = corresponds to the bosonic case, q = 1 to the 
fermionic case), and Ai the statistical vector potential. An harmonic potential term has 
been added as a long distance regulator. It discretizes the energy spectrum, and allows 
for a convenient computation of thermodynamical quantities in the thermodynamic limit 
(w ^ 0, where lo is the harmonic oscillator pulsation). 

In a recent series of papers [3], the A^-anyon partition function has been shown to be 
written as 

Z^(a)=Ei^H«)n^(^)" (2) 

where summation is understood over all partitions V of N . Each partition ^^^ Lvi^ = N 
is in one-to-one correspondence with an equivalence class of the permutation group S'at, 
namely the class of permutations which can be written as a product of cycles of length L, 
V]^ times. The one-particle partition function Zi{l3) = l/(2sinh|)^ is obviously statistics 
independant { ^ = h/3uj,P = ^). Note that in (2), the statistical information has been 
completely factorized out in the F-p(a)'s, which satisfy 

F-p{l + a) = sgn{V)F-p{a) (3) 

It follows that Zp4{a) is naturally written as a sum of symmetric and antisymmetric terms 
under a — > 1 + q. By convention one has F-p(O) = 1 for all V, i.e. q = corresponds to 
Bose statistics. 



In a brownian motion language, the F-p's are simply expressed as 

Fp(a)=<e*"S/<^^" >(C) (4) 

where the average < >fc) i^ taken over all possible brownian paths which induce, in the 
configuration space, a given permutation V, i.e. which correspond to a given partition 
of the integer N. Note that these curves are closed curves in the configuration space 
quotiented by the permutation group Sn- The 0/j's are the relative angles between any 
pair of particles {I, J} (see below for a more precise definition). For example, F/v is 
associated with cycles of maximal length L = N, with i^l = 1. Such cycles, involving the 
N particles of the system, are by definition A^-particles paths over a time interval r = hp. 
Equivalently, they can be considered as 1-particle paths over a time interval Nt, simply 
because the particles are identical [3]. Therefore, < >(c) can be taken over all brownian 
paths of length Nt. These paths are divided into N segments of equal length r, each 
segment corresponding to the motion of a point Mj (all the M/'s have the same speed). 
The 0/j's are defined as the angles of the vectors MjMj with a fixed axis in the plane 
( — oo < 9jj < +oo) (see fig. 1 for F3). 

In the simplest N = 2 case, eq. (2), has been exactly solved [1]. It is instructive to 
display the results and their intimate connection with probability distributions of planar 
brownian windings. .^2(0) is given as 

Z2{a) = ^F2(a)Zi(2/3) + ^Fn{a)Zi{pf (5) 

and is computed from the complete linear states 

V'n,™(r, e) = e^-^^e-Z^^^I— "l^™-"! (^^2^ (6) 

En,m = i-^C^n + \m - a\ + I) (7) 



where the L^'s are the Laguerre polynomials. One gets for F2{a) and -Fii(a) 



sinh(l — 2a)| cosh(l — 2a)| 

sinh 4 cosh 4 



^2(«) = \ . ^ '' Fnia) = ^ . '' (8) 

I cosh J 



In terms of brownian winding, the Fourier transform of Fii{a) (0 < a < 1) yields the prob- 
ability for a brownian curve to wind an even number of half-windings around a fixed point, 
i.e. the usual result for an integer brownian winding [4]. ^2(0) leads to the probability 
for an odd number of half-windings [3]. 

The second virial coefficient follows directly [5] 

1 {1-af 
a2[a) = - (9) 

One notes that in the thermodynamic limit where the harmonic potential vanishes 

F2ia) ^^^0 F^^\a) = 1 - 2a (10) 



Fn{a) ^^^0 Fi;^\a) = 1 (11) 

The trivial thermodynamical limit of Fii{a) is easily understandable in terms of (4), since 
in this limit there is no chance for the 2 independent 1-particle brownian paths to overlap. 
In the A^ = 3 case, things become non trivial, since only a part of the spectrum (the 
linear states [6] that generalize (6)) is known. One has by definition 

Z3ia) = ^F3(a)Zi(3/3) + iF2i(a)Zi(/3)Zi(2/3) + ^Fiii(a)Zi(/3)3 (12) 

6 2 b 

It happens that the antisymmetric part of .^3(0), namely F2i{a), is known [7] to be exactly 
built from the linear states, i.e. 

Zi(/?)Zi(2/3)F2i(a) = Zf"(a) - Zf"(a + 1) (13) 



with [6] 



This yields 



with 



7/m/ N cosh[3(l - a)C] 



^'^^""^ 32sinh2|sinh2Csinh2f ^^"^^ 



, , sinh(l-2a)f , , 

^21 (a) = -rr^ (15) 

smh -^ 



F2i{a) ^^^0 Fi^\a) (16) 



as it should. 

In [3], numerical estimations for F^ (o), and for Fiii{a) at order ^^, have been per- 
formed in the thermodynamic limit 

Ff(a)^(l-3«)(l-^«) (17) 

This numerical result, altogether with F^^((q), is sufficient to propose a simple analytical 
expression for the third virial coefficient [3] 

as (a) ~ ^ + Y^ sin2(7ra) (18) 

Having reached this point, and looking at the simple expression for F^ (ct), it seems 
natural to suspect that it could be generalized to iV > 3, and that it might be derived 
from pure theoretic considerations. This is what is now going to be shown. 

Indeed, using (4), one is tempted to make numeric simulations on a lattice for the 
F-p(a)'s. Note that in the thermodynamic limit, one has always that the Fj, (a)'s factor 
out in a product of Ff^ (a)'s, where the L's enter in the partition V. From this consider- 
ation, it follows that only the Fpf{a)^s should be easily accessible numerically. Note also 
that the regularisation due to the lattice spacing necessarily alters the results, since two 
points cannot approach each other arbitrarily close ([8,9]). In the case of a true brownian 



motion, when two points Mj and Mj come close to each other, the accumulated angle Ojj 
has a non negligible probability to become arbitrarily large (the typical brownian path 
being non-differentiable, 6jj follows a broad law probability distribution). Thus, when 
one tries to simulate brownian paths on a lattice, a correction is necessary to take into 
account this characteristics of brownian winding. Each time two points Mj and Mj ran- 
domly walking on the lattice are separated by only one lattice spacing, one adds to Ojj a 
multiple of 2tt which is drawn accordingly to a Cauchy law [9]. Typically, we have used 
random walks ranging from 50000 to 150000 steps. The resulting numeric simulations for 
the -F^ (a)'s suggest the following simple formula 

N-l ^ 

Fj^\a)=l[il-ja) (19) 

fc=i 

valid for < a < 1 (to obtain an expression in the whole interval [0, 2] one uses F^ (a) = 
FP{-a) = Fi^\a + 2)). 

Eq.(19) correctly reproduces the N = 2 and A^ = 3 cases ; it also has the right a = 
and a = 1 limits ; finally it satisfies -FAr(l + a) = {—l)^~^^F]\^{a). Note that a similar 
formula has recently appeared in an a priori different context [10], the one-dimensional 
Calogero model [11]. There, the -F^ (a)'s, are derived from the exact one dimensionnal 
A^-body spectrum (they are the coefficients appearing in the expansion of the density in 
powers of the fugacity). This similarity is not a surprise : it is well-known that there is a 
correspondence [12] between the Calogero model and the anyon model projected into the 
lowest Landau level of an external magnetic field. This might imply that when anyons are 
projected into the Landau groundstate, mostly cycles of maximal length contribute to the 
thermodynamics. 

The expression obtained for F^ (a) is again quite simple. One should be able to find 
a way to support it by a pure analytical approach. Let us concentrate on the A^ = 3 case. 



As already stated, the antisymmetric part F2i{a) of Z^i^a) is determined by the known 
hnear states ; its symmetric part is 

Zs{a) + Zs{l + a) = ^F3(a)Zi(3/3) + ^Fin{a)Zi{pf (20) 

Could it be that -F3 is also determined by the known part of the 3-anyon spectrum ? Some 
simple considerations [13] on Fm will show us the way. The paths contributing to Fiii(a) 
are simple one-particle cycles going back to their initial position at the time r = hp. The 
final configuration is identical to the initial one, thus there is no need for quotienting the 
configuration space by the permutation group 53. It follows that Zi{(3)^ ; Fiii(a) is simply 
the partition function of a system of non-identical particles (Boltzmann statistics) inter- 
acting topologically a la Aharonov-Bohm with coupling constant a now defined modulo 1 
instead of modulo 2 in the identical particle case. In the case of 2-anyon one indeed has 
that 

Zi{(3f ■,Fnia) = Z2{a) + Z2{a + l) (21) 

meaning that Fii{a) (for 2 non identical particles) can be directly obtained from ^2(0) 
(the partition function for 2 bosonic identical particles)^ Thus one writes for Fm 

Zi{(3f ;Fiii(a) = ^3(0, a, a) +^3(0 + l,a + l,a + l) +2Z3{mixed) (22) 

with 

Z^^mixed) = Zj,(a, a, a + 1) + ^3(0;, a + 1, a + 1) (23) 

where one has explicitly denoted the relative statistical coefficient for each pair of particles 

(here a or a + 1). In the case a = 0, this implies that the mixed representation of 5*3 is 
■^In other words in ^2(0) one has summed on even angular momentum quantum numbers, and, in order 
to get the complete spectrum, one simply shifts a — > a + 1 to get the odd angular momentum quantum 
numbers. 



represented either by (0, 0, 1) or by (0, 1, 1), meaning in the first case that the pairs {1, 2} 
and {1, 3} are bosonic, and {2, 3} is fermionic, and in the second case that the pair {1, 2} 
is bosonic, whereas {1,3} and {2,3} are fermionic. The factor 2 in front of Z^{mixed) 
is due to the fact that the symmetric and antisymmetric representations of S^ are one- 
dimensional, but the mixed one is two-dimensional implying doubly degenerated states in 
the spectrum. 

This in turn implies for F2,{a) 

Zi (3/3)^3(0) = ^3(0, a, a) + Z^{a + 1, a + 1, a + 1) — Z^{m,ixed) (24) 

This formula is analogous to (13), so one might ask about its validity when one restricts 
oneself to linear states. One has the generalized linear eigenstates [14] 

KJ 

-E„,r„,j=^(2n+^|m/j-a/j| + iV-l) (26) 

KJ 
for a topological Hamiltonian where the Aharonov-Bohm coupling constants ajj depend 

on the pair {/, J} of particles coupled {p^ = J2i<j''''jj)- Mixed linear states are obtained 

by trading for each couple {/, J} of particles the statistical parameter a for a statistical 

parameter ajj. One can easily show that the partition function for such mixed linear 

states is then directly obtained from the linear partition function Z3 (a) by — ^"-2 — -a ; — > 

Y^KjCtiJ- For iV = 3 

Z3 (a,a,l + aj " 32 sinh^ | sinh^ g sinh^ f .27) 

4-(a,l + a,l + a) = '^osh[(i 3a)g] 

■^ ^ ' ' ^ 32 sinh2 I sinh2 § sinh^ ^ 

Let us now consider (24) restricted to linear states only. One finds that if the expression 
Z^*"(q, a, a) + Zf "(a + 1, a + 1, a + 1) - Z^*"(a, a, 1 + a) - 4*"(a, 1 + a, 1 + a) (28) 



has not the correct finite limit -F3 (a) in the thermodynamic limit, but on the contrary 
diverges, the left over divergence does not depend on a and simply reads Zi{(3)Zi{2>(3). It 
is quite a striking fact that when one adds Zi(/3)Zi(3/3) to Z-^"^{a, a, a + 1) + Z^"'{a, a + 
1, a + 1) in (28) one gets in the thermodynamic limit exactly the right F^ (a) as given in 
(17). Moreover, one can compute F-i{a) as given by (28) for an arbitrary ^ 

F, («) = ^mh(2-3a)|sinh(l-3a)| ^^^^ 

sinh I sinh ^ 

This expression generalizes quite nicely F2{a) given in (8). 

Similar considerations (albeit less transparent) can be made for the Fjy s, A^ > 3. 
One notes that a possible generalization of (29) 

^ ^ ^Vsinh(A;-iVa)|\ , ^ 

would automatically satisfy the right thermodynamic limit (19) obtained in the numerical 
simulations. 

To conclude this analysis, let us come back to the standard case a = 0. One knows 
that Zf{(3)Fiii{0) = Zf{P) simply describes 3 independent bidimensionnal oscillators ; 
one also has 

.^rnnm- cosh(30 + 2cosh^(g/2) 

^^^°' °' "^ - 32sinhn,inh2esinh^M ^''^ 

l + 2cosh2(^/2) 
^'^^'^'^^"32sinh2|sinh2esinh2f ^^'^' 

One would like to find simple expressions for ^3(0,0, 1) and .^3(0, 1, 1). It is clear from 

the bosonic (31) and fermionic (32) 3-body partition functions that the linear part of 

the spectrum contributes respectively to cosh(3^)/(i and 1/d (i.e. cosh(3 — 3a)S,/d for 

a = 0, 1) whereas the common 2 cosh {^/2)/d is built by the unknown part of the spectrum 



((i = 32sinh | sinh ^sinh ^). It is quite natural to try for ^3(0,0, 1) and ^3(0, 1, 1) 

zmnn cosh(20+2cosh^(g/2) 
^^^°'°''^'~^32sinh^|sinh^esinh^f ^''^ 

7 m 1 n cosh^ + 2cosh^(^/2) 

^^^°''''^'^32sinh^|sinh^esinh^f ^''^ 

and to check if (22) holds for a = 0. One finds again the same left over divergence 
Zi{l5)Zi{?>l3) as above (eventhough one gets the correct behavior --^ i in the thermody- 
namic limit), that should be added to ^3(0, 0,1) + ^3(0, 1,1) to get Zi{f3)^. It would 
certainly be interesting to understand more precisely, already at the level of (22), the 
origin of this additional term, and also in which way it splits between Z3 (0,0,1) and 
^3(0, 1, 1), thus allowing a complete knowledge of the 3-body partition functions for the 
mixed eigenstates of S3. 

Acknowledgements : J. D. and S. O. acknowledge stimulating conversations with J. 
Myrheim, and also for drawing our attention on [10]. 

Figure Captions : 

A brownian closed path contributing to F^ («)• The points A,B and C divide the 
curve into 3 patches of equal length. The 3 points Mi, M2 and M3 span each of the 
patches with an equal speed (Mi goes from A to B, etc.). 
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